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Abstract. This article reviews on-shell methods for analytic computation of loop 
amplitudes, emphasizing techniques based on unitarity cuts. Unitarity techniques 
are formulated generally but have been especially useful for calculating one-loop 
amplitudes in massless theories such as Yang-Mills theory, QCD, and QED. 
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1. Introduction 

Scattering amplitudes can be constructed in terms of their singularities. For tree 
amplitudes, these are complex poles. In loop amplitudes, there are branch cuts, as well 
as other singularities associated with "generahzed" cuts. All of these singularities probe 
factorization limits of the amplitude: they select kinematics where some propagators are 
put on shell. Thus, the calculation can be packaged in terms of lower-order amplitudes 
instead of the complete sum of Feynman diagrams. 

The "unitarity method" started as a framework for one-loop calculations. Instead 
of the explicit set of loop Feynman diagrams, the basic reference point is the linear 
expansion of the amplitude function in a basis of "master integrals," multiplied by 
coefficients that are rational functions of the kinematic variables. The point is that the 
most difficult part of the calculation, namely integration over the loop momentum, can 
be done once and for all, with explicit evaluations of the master integrals. The master 
integrals contain all the logarithmic functions. It then remains to find their coefficients. 

If an amplitude is uniquely determined by its branch cuts, it is said to 
be cut-constructible. All one-loop amplitudes arc cut-constructible in dimensional 
regularization, provided that the full dimensional dependence is kept in evaluating the 
branch cut. Each master integral has a distinct branch cut, uniquely identified by its 
logarithmic arguments. Therefore, the decomposition in master integrals can be used 
to solve for their coefficients separately using analytic properties. It is not necessary to 
reconstruct the amplitude from the cut in a traditional way from a dispersion integral. 
Rather, we overlay information from various cuts separately. 

Beyond one-loop order, not many analytic results have been obtained by unitarity 
methods in the absence of maximal supersymmetry. Perhaps the main obstacle is the 
lack of a manageable, systematic, explicitly known set of master integrals. Also, D- 
dimensional unitarity cuts of higher-loop amplitudes involve lower-order amplitudes 
which still contain loops and yet have D-dimensional momenta on some external legs 
(the cut lines). 

Analytic calculations are simplest in massless theories, where formulas can be 
written compactly in the spinor-helicity formalism. Spinor variables are helpful inside 
the loop as well, when propagators associated to massless field are placed on shell in a 
unitarity cut. For this reason as well, we work mostly in four-dimensional Minkowski 
space and its analytic continuations. 

This article is organized as follows. In Section 2, we introduce the master integrals 
and review integral reduction. In Section 3, we describe the unitarity method and cut 
integrals. We evaluate the cuts of master integrals, explain the evaluation of general cut 
amplitudes, and list formulas for the coefficients of of master integrals, given a general 
one-loop integrand. In Section 4, we discuss generalized unitarity cuts for one-loop 
amplitudes, from quadruple and triple cuts to single cuts along with adaptations of 
Feynman's Tree Theorem. In Section 5, we address D-dimensional unitarity methods, 
which can be used either for an exact expression of the regulated amplitude or to 
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Figure 1. One- loop master integrals: box, triangle, bubble and tadpole. The Lorentz 
vectors are sums of external momenta, all directed outward. 

avoid the separate calculation of "rational terms" in addition to the master integral 
coefficients. In Section 6, we mention issues that arise in the presence of massive 
particles. In Section 7, we review applications of MHV-diagram constructions. In 
Section 8, we discuss on-shell recursion relations for coefficients or rational terms in 
one- loop amplitudes. In section 9, we summarize recently published analytic results for 
one- loop amplitudes. Finally, in section 10, we look at prospects for continuing progress 
along these lines to higher-loop amplitudes. 



2. One-Loop Master Integrals 

One-loop Feynman integrals can be expressed as a linear combination of master integrals 
with rational coefficients. The master integrals are typically taken to be the set of scalar 
integrals, meaning that there is no tensor structure left in the numerator. The scalar 
n-point integral is defined as 

d'^i 1 



/„ = (-ir+iz(47r)f J 



(27r)^ (£2 - ml){{i - K^f - ml){{i - - K^f - ml) ■ ■ ■ ((£ + KnY - ml] 
Here, the Ki are sums of external momenta. They are strictly four- dimensional. 

Working in theories with mostly massless particles, it is natural to use dimensional 
regularization, i.e. Z) = 4 — 2e. We are often satisfied with a result good to C(e). In 
analytic computations, we keep the external momenta fixed in four dimensions in order 
to use the spinor-helicity formalism. 

Integral reduction [Il|2l|3llH[5l|6l[7] is a clearly defined procedure for expressing 
any one-loop Feynman integral as a linear combination of scalar boxes, scalar triangles, 
scalar bubbles, and scalar tadpoles, with rational coefficients: 

^l-loop^^ E Cn(K)/n(K) (1) 

n K={/(-i,...,if„} 

In four dimensions, n ranges from 1 to 4. In dimensional regularization, the tadpole 
contributions with n = 1 arise only with internal masses. If we keep higher order 
contributions in e, we find that the pentagons [n = 5) are independent as well. 

Let us briefly outline the traditional reduction procedure. We assume that the 
integral has been constructed from Feynman diagrams, so that the denominators are 
propagators of the form Di = {i — Pi)^ — Mf, along with the propagator defining the 
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loop momentum, Dq = P — Mq, chosen anew for each term at each stage. There are 
three steps. 

First, we ehminate tensor structure (i.e. momentum-dependent numerators) in 
terms with at least five propagators. Any appearance of £^ in the numerator is replaced 
by Mq + Do, and the Do term cancels against the denominator. The remaining 
momentum dependence in the numerator is polynomial in contractions of the form i ■ P. 
Among the five propagators, there are four independent momentum vectors Pi in which 
to expand any P. Then we make the replacement 2i- Pi = + P^ — Mf + Dq — Di and 
cancel Do and Di against the denominator. Step by step, the degree of the polynomial is 
lowered until we have a scalar numerator or at most four propagators in the denominator. 

Second, we eliminate remaining tensor structure in the terms with at most four 
propagators. This is done by using the momenta appearing in the denominators to build 
a basis of Lorentz-covariant tensors in which to expand the integral. Contracting the 
tensors with external momenta gives the constraints needed to solve the linear system. 
It can be particularly efficient to use contractions with complex momenta constructed 
from spinors associated to different external legs [H [9], [TOl IHl [12] . 

Third, we express n-point scalar integrals with n > 4 in terms of lower-point scalar 
integrals. If n > 6, then there is a nontrivial solution {oj} to the five equations 



Er=i«i = and YJl=i(^iPt = 0. With this solution, «i A = E^«i(^' - Mf). 



Divide the integrand by the (momentum-independent) right-hand side of this equation 
and multiply it by the left-hand side. The factors Di will cancel against the denominator 
and reduce n by one. The final remaining concern is the scalar pentagon. If we are 
keeping full e dependence in dimensional regularization, it is an independent master 
integral. If we truncate the integrals at 0(e), then the scalar pentagon can be reduced 
further to four scalar boxes. The scalar pentagon integral is finite, so we can now 
treat its loop momentum as four-dimensional. The final reduction involves expanding 
it in terms of the axial vectors constructed from triples of four independent external 



Explicit formulas for the divergent one- loop master integrals are given in [13], for 
all arrangements of massive and massless lines. In the next few sections, we will study 
unitarity cuts specifically at the case in which all internal lines are massless. Then the 
tadpoles are absent. We now list the necessary master integrals in this case evaluated 
through 0{e). The expressions here are taken from [TU[T5]. Other useful expressions for 
scalar box integrals, convenient for analytic continuation to different kinematic regions, 
appear in [T6] . 

The dimensional regularization parameter is e = (4 — D)/2. The constant rr is 
defined by 



momenta. 




(2) 



Scalar bubble integral, no internal masses: 




(3) 
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Scalar triangle integrals, no internal masses: 

II Kl — Kl — {) and Kl^Q, then the scalar triangle is called "one- mass", and it is 



jlm 
-'3 



If i^l = and Kl, K2 7^ 0, then the scalar triangle is called "two-mass", and it is 



j2m 

-'3 — 



rr {-Kl 



The "three-mass" scalar triangle is finite and given by 

3 



T3m 
-'3 



E 



Li-; 



i6j 



1 — i6-j 



— Lio 



t6j 



i6i 



+ 0{e), 



where we have defined the following: 

A3 = -(Kf)^ - {Klf - {Klf + 2KlKl + 2KlKl + 2KIkI 
^ _2K]-{K! + Kl + Ki) 



(4) 
(5) 

(6) 

(7) 
(8) 



Scalar box integrals, no internal masses: 

Let s = {Kl + K2)^ and t = {Ki + K^y. The dilogarithm function is defined by 
Li2(a;) = — /g^ln(l — z)dz/z. 

If all four momenta are massless, i.e. Kf = Kl = K^ = Kl = (a special case for 
four-point amplitudes), then the box integral is given by 



rOm 
-'4 



rr 
st 



(9) 



If only one of the four momenta, say Ki, is massive, and the other are massless, i.e. 
Kl — K^ — Kl — 0, then the box is called "one-mass" , and it is given by 



rim 
-'4 



2rr 1 

~^ ^ 
2rr 
~ It 



Li2 1 - 



+ Lis 1 



t 



iln^ f£ 
2 \t 



+ 



~6 



(10) 
+ 0{e). 



There are two distinct arrangements of two massive and two massless legs on the corners 
of a box. In the "two-mass-easy" box, the massless legs are diagonally opposite. If 
Kl = K'l = while the other two legs are massive, the integral is 

2rr 1 



T2m e 

-'4 — 



st 



2rr 



-Kir] 



11) 



st - KlKl 



Li2(l-^)+Li2 (l-^ 



+ Lis 1 - 



Kl 
t 



- Li, 1 - 



KlKl 



st 



t 



-In^ 
2 \t 



Li2 ( 1 - ^ 



+ 0{e) 



In the "two-mass-hard" box, the massless legs are adjacent, li Kl — Kl — Q while the 
other two legs are massive, the integral is 



r2m h 



2rr 1 



l{-s)-^ + {-t)-^-l{-Klr-\{-Kl) 



2\-e 



(12) 
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2rr 
It 



If exactly one leg is massless, say = 0, then we have the "three-mass" box, given by 



2rr 



st-KfKi e2 
2rr 



(13) 



Finally, the "four-mass" box, which is finite, is given by 

2 



4m 

St 



jAm 
-'4 



1 



- L12 1 + ^— 

— S — 26 

—t — is 



Xj -7] 



^ Mn IH ^ — — : 



- Li2 1 



+ \n{-Xj){\n{-Kl - le) + ln(-s - le) - \a{-Kl - le) - \n{-Kl - le))) . 

Here we have defined 

T]{x,y) = 2TTi[i){—lm Im y)i){lm (xy)) — 'i9(Im x)'i9(Im y)'d{—lm (xy))], 

and Xi and X2 are the roots of a quadratic polynomial: 

ax^ + bx + c + ied = a{x — xi){x — X2), (14) 

with 



— s — te 
—t — is 



In 1 



-s — le 
-t — is 



a = tKi, b = st + KlKl - K^Kl c = sKl 



d 



-Kl 



2.1. Simple theories 

Some field theories are known to obey additional constraints on their expansion in master 
integrals. Maximally supersymmetric theories, Af = 4 SYM and J\f = 8 supergravity, 
have a "no-triangle" property [T71[T8],|19l[20l[2T],[22]. One-loop amplitudes can be written 
in terms of box integrals only, without triangles or bubbles. This property follows from 
supersymmetric cancellations combined with power-counting, although for supergravity, 
new reduction formulas needed to be found. Master integrals for higher-loop amplitudes 
should not contain triangle or bubble subdiagrams, either. 

Massless supersymmetric theories should be four-dimensionally cut-constructible. 
That is, the expansion in master integrals ([H) is valid when truncated at C(e), without 
any additional rational terms that could arise from the expansion in e. This was shown 
specifically for color-ordered gauge theory amplitudes [171 [H] but is expected to be 
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Condition (C): TrR(n7^iT"') = m Tradj(n7=iT"»), n<p 


Non-susy theories have 


only boxes 


no bubbles 


a Rf satisfies C with 


p=6, m=4 


p=4,m=4 


and Rs satisfies C with 


p=6, m=6 


p=4,m=6. 


Susy theories have 


only boxes 


no bubbles 


if R^ satisfies C with 


p=5, m=3 


p=2,m=3. 



Table 1. Conditions for cancellations of bubbles and triangles. Matter representations 
are denoted by / for fermions, s for scalar and x for chiral supermultiplets. Taken from 

more generally valid, and is also useful in organizing calculations in nonsupersymmetric 
theories. 

Certain gauge theories with reduced supersymmetry still feature enough symmetry 
to guarantee the cancellation of triangle or bubble contributions. In A/" = 6 supergravity, 
bubble contributions are absent at one loop [23]. In Yang-Mills, theories with triangle 
or bubble cancellations have been recently characterized in terms of the representation 
of matter content in |23] , whose conclusions we list here in Table 12.11 An example of 
a theory with box integrals only is A/" = 2 with one hypermultiplet transforming in 
the symmetric tensor representation of the gauge group SU{N) (where > 3) and 
another hypermultiplet transforming in the antisymmetric tensor representation. An 
example of a theory without bubbles is A/" = 2 with an SU{N) gauge group and 2A^ 
hyp ermult iplet s . 

Multi-photon amplitudes in QED at one loop exhibit simplicity as well. Furry's 
theorem implies that the n-photon amplitude vanishes if n is odd. Regarding the integral 
basis, it has been shown [25] that bubble integrals are absent for n > A, and both bubble 
and triangle integrals are absent for n > 6. Rational terms appear only for n = 4. 

3. Unitarity methods 

The unitarity cut of a one-loop amplitude is its discontinuity across the branch cut in a 
kinematic region associated to a particular momentum channel. The name comes from 
the unitarity of the S-matrix: since S'^S = 1, and we expand S = 1 + iT where T is the 
interaction matrix, then 2Im T = T^T. Expanding this equation perturbatively in the 
coupling constant, we see that the imaginary part of the one- loop amplitude is related 
to a product of two tree-level amplitudes. Effectively, two propagators within the loop 
are restricted to their mass shells. This imaginary part should be viewed more generally 
as a discontinuity across a branch cut singularity of the amplitude — in a kinematic 
configuration where one kinematic invariant, say K"^, is positive, while all others are 
negative. This condition isolates the momentum channel K of interest; K is the sum 
of some of the external momenta. We will take cuts in various momentum channels to 
construct the amplitude. 
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Figure 2. Unitarity cut of a one-loop amplitude in the K momentum channel. The 
two propagators are constrained to their respective mass shells. The disks represent 
the sum of all Feynman diagrams linking the fixed external lines and the two cut 
propagators. 



For a one-loop amplitude, the value of the unitarity cut is given by Cutkosky rules 
[26]l|l The Cutkosky rules are expressed in the cut integral, 

A^'"^°°^ ^jdfi A'-J, X „ (15) 
where the Lorentz-invariant phase space (LIPS) measure is defined by 

dfi = d% d% 5(^)(^i +i2-K) 6^+\il) 6^^\£l). (16) 

Here, the superscript (+) on the delta functions for the cut propagators denotes the 
choice of a positive-energy solution. 

How can these unitarity cuts be used to calculate the amplitude? By applying 
the cut A in various momentum channels, we get information about the coefficients of 
master integrals. 

Consider applying a unitarity cut to the expansion ([T]) of an amplitude in master 
integrals. Since the coefficients are rational functions, the branch cuts are located only 
in the master integrals. Thus we find that 

AA'-ioop^J2 c„(K)A4(K) (17) 

n K={Ki,...,K„} 

Equation (fT7|) is the key to the unitarity method. It has two important features. 
First, we see from (IT^ that it is a relation involving tree-level quantities. Second, many 
of the terms on the right-hand side vanish, because only a subset of master integrals 
have a cut involving the given momentum K. Meanwhile, we enjoy the freedom of using 
all possible values of K in turn. In effect, we have traded the original single equation 
([1]) for a system of several shorter equations. 

Our remaining task is to isolate the individual coefficients Cj. With generalized 
unitarity, quadruple cuts isolate the box coefficients. We have to work a little harder to 

I This technique and related approaches are discussed in |27| in their original context, which excluded 
massless theories. Modern interpretations were introduced in [iTl fT4]. 
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get the triangle and bubble coefficients. In the following section we will discuss triple 
cuts for triangle coefficients. Here, we recall that four-dimensional cut-constructibility 
allows us to distinguish box, triangle, and bubble contributions from standard unitarity 
cuts, and we will see how these coefficients can be calculated in practice. 



3.1. Cuts of master integrals 



The utility of equation (jlSil depends on knowing the master integrals and hence their 
branch cuts. Having listed the master integrals for massless theories in the previous 
section, we can calculate the cuts explicitly by taking the imaginary parts of these 
functions in various kinematic regions. Recall that the kinematic region associated to 
the unitarity cut in momentum channel K has > and all other invariants negative. 
In the limit of large K"^, we see that the behavior of the master integrals features uniquely 
identifiable products of logarithms [13]. It follows that no linear combination of master 
integrals with rational coefficients, in a given momentum channel, can be cut-free. 

The cuts of the bubble integrals are purely rational; this is easily seen from the 
ordinary logarithm in equation (jS]). The cuts of all other master integrals are logarithmic. 
The various arguments of the logarithms identify the original master integrals. 

We need to evaluate the left hand side of equation ( fTTl) . by carrying out the 2- 
dimensional integral of ( lT5l) . This is neatly accomplished by the Cauchy residue theorem, 
in a technique known as "spinor integration" [281 129] . We now illustrate the technique 
applied to the master integrals themselves, starting from their definition in ([T]). 

To implement the cut conditions, it is convenient to reparametrize the loop 
momentum in terms of spinor variables. Now it is crucial that our cut is in 4 dimensions. 
(The D-dimensional generalization will be presented in Section |5l) 

Since ii is null, we can parametrize it with 

{il)aa = tXa~Xa, (18) 

where A^, are homogeneous spinors (taking values in CP^), and t takes nonnegative 
values. The original loop momentum is real-valued, so we will integrate over the 
contour where A is the complex conjugate of A . In the integral measure, we make 
the replacement [30 j§l 



A dX 



(19) 



d^ii 5^+) (£')(•) = - / -dt (A dX) 
Jo 4 Jx=x 

Now we make this substitution explicitly in the second delta function of the LIPS 
measure defined in ( flGll . The momentum of the second cut propagator is {£2)ad = 



tXaXa, SO the measure becomes 



Idt 
4 



dfi 



(A dX) 



x=x 



X dX 



tiX\K\X 



)(.). (20) 



§ The factor of —4 has been omitted in much of the hterature, but such constant factors cancel out in 
unitarity methods and can be disregarded. 
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This second delta function sets t to the value 

t -- 



X\K\\ 



so, taking account of the prefactor (Alf^lA 
perform the t-integral trivially: 



(A d\) 



(21) 

of t inside the delta function, we can now 



A dX 



A(X\K\X 



(22) 



The remaining integral over the spinor variables is carried out with the residue theorem. 
We will see how this is done in the master integrals before proceeding to the case of 
general amplitudes. 

3.1.1. Cut bubble Let us start with the scalar bubble integral. The integrand consists 
entirely of the two cut propagators, so the cut is simply the integral of the LIPS measure. 



A 



1 



£2(£ - Kf 



dfi 



(A dX) 



x=x 



X dX 



4(A|K|Ar 



(23) 



In calculating coefficients of complete amplitudes, it can suffice to leave the cut bubble 
in the form (l23l) and work at the integrand level. Here we continue and show how 
to apply the residue theorem to complete the integral. We make use of the following 
identity to rewrite the integrand as a total derivative [30]. Here r] is an arbitrary spinor. 



[A dX] 



X\K\X 



[dX %] 



[A V] 



{X\K\7^](X\K\X 



(24) 



However, the integral is not identically zero, because there are delta-function 
contributions along the contour. In the theory of a complex variable, we know that 
d 1 



27r6{z-b). 



(25) 



dz {z — b) 

Therefore, we pick up a residue at the pole |A) = l-ft'lr]]. Along the contour, since A and 
A are conjugates, we also substitute |A] = \K\ri). The result for the four- dimensional 
cut bubble is thus 



Ah 



:A 



' 27r 



[A V] 



X\K\X 



1 

27(1' 



(26) 



|A> = |i^|r?] 

(Different conventions in the literature such as in [28] and [31] yield results with different 
powers of i and 27r; again these will be unimportant as long as the framework is 
consistent. The main point is clear, namely that the result is a constant, proportional 
to the discontinuity across the branch cut of the logarithmic function dH]). ) 
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3.1.2. Cut triangle In the unitarity cut of the scalar triangle, there is one propagator 
left over along with the LIPS measure. Converting to the spinor variables, this factor is 

+K'l- Performing the t integral as before, making the substitution 



{l+K^Y = t{\\K^\\ 
(!2T|) throughout, we have 
1 



A 



where 



(A dX) 



A dX 



4(X\K\X 



X\Q\X 



(27) 



(28) 



Again, it is worth leaving the expression in the form (^7^, but let us see how to finish 
the integral. The two factors in the denominator can be combined with a Feynman 
parameter, and the spinor integral done just as in the bubble case, so that we have 



dx 







(A dX) 



A=A 



A dX 



A{X\{l-x)K + xQ\X 



x)K + xQy 



.(29) 



The result for the cut in the i^-channel is 

1 



Ah 



4a 

7r^ 



In 



-2iKi + K ■ K,) + ./^s 



where A3 is defined in (^^. Notice that this result is logarithmic, as expected. Moreover, 
it is clear that all three-mass triangles are uniquely identified by the functions A3, which 
play a distinguished role in the expression as the arguments of square roots. (For one- 
mass and two-mass triangles, A3 is a perfect square, so the square roots disappear from 
the formula while the logarithm remains.) 



3.1.3. Cut box The calculation of the cut scalar box integral is similar. Now there are 
two uncut propagators identifying the box, which we write as {£ — and (£ — Kj^. 
Converting to the spinor variables, they become K"^ — t (^X\Ki\X and K'j — t(^X\Kj\X , 
respectively. Performing the t integral and making the substitution fl2T|) throughout, we 
have 

1 



A 



1 



(A dX) 



X dX 



^ 4K^X\Q,\X (X\Q,\X 



(31) 



where now we define Qi,Qj by 



(32) 



Here again, we can evaluate the integral by introducing a Feynman parameter. It takes 
a form similar to the triangle. The final result is 

'Qi-Qj + ^fKi 



Ah 



:ln 



/A 



(33) 
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where 

A,, ^ (Q,, ■ Q^y - QlQ]. (34) 

We see that the cut is again logarithmic. One can check that for either of the two choices 
of the cut configuration (straight across the box or selecting one corner), the function 
Ajj under the square root corresponds to the discriminant of the quadratic polynomial 
in (JT^ . In the cases where any of the corners of the box is a null momentum, A^- is a 
perfect square. 



3.2. Unitarity cut of the amplitude 

To evaluate the cuts of master integrals, we needed only the relatively simple identity 
(^^, from which we could find the residues according to (!25|) . The cut of the full 
amplitude will generally have more complicated dependence on the loop momentum. 
Let us examine its analytic structure. 



3.2.1. Stokes' Theorem on the complex plane Viewed formally, the residue theorem 
needed to evaluate cut integrals is an application of Stokes' Theorem to the complex 
plane [31]. Indeed, it is intuitively helpful to rewrite the cut integral in the familiar 
language of a single complex variable. Given the cut momentum K, we can choose two 
null momenta p and q such that 

K = p + q. (35) 

To see this, choose an arbitrary null vector p, and solve the equation q = K — op such 
that q^ = 0. The solution is a = {p ■ K)/{2K'^). Then set p = ap. 

The spinors A, A parametrizing loop momentum can now be expanded in the basis 
of spinors from p and q 

A = Ap + zXq, A = Ap + zXq. (36) 

The homogeneity of the spinors has been used to set the coefficients of Ap and Ap to 1, 
and now we see the familiar complex variables z and z emerge as the other coefficients. 
Recall that A, A were coordinates on CP^. Equation (!36|) is a representation of the 
standard mapping of the complex plane to a sphere. 

The cut integral measure now takes a reasonably simple form, since (Ac/A) [Ac/A] = 
—K'^dzdz, and ( AI-ft^lA 



K^{1 + zz). Thus the variable t gets fixed to the value 

t = (37) 
1 + zz 



and the cut measure fl22l) becomes 



The cut integration can be performed by the Generalized Cauchy Formula. Given 
an integrand F{z,z), construct a primitive G{z,z) with respect to z. Let D be a disk 
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in the complex plane containing all poles in z of G{z, z). Then 

/ F{z,z) dz Adz = f dzG{z,z)-2m V Res{G{z,z),Zj}. (39) 

Jd JdD . . , 



poles Zj 



The bubble contribution is distinguished by being rational. 



3.2.2. Partial fraction identities for splitting denominator factors To look a little 
deeper into the structure of cuts of amplitudes, we return to the spinor variables A, A. 
Suppose we take the expressions for the tree amplitudes in the cut integral ( fTSll from the 
Feynman rules. The integrand is a rational function whose denominator is a product of 
propagator factors of the form (£ — KiY. As we have seen the cuts of master integrals, 
such a factor becomes 



Other factors of ( Ali^lA 



X\K\X 



(40) 



arise from the integral measure and the substitution for t 
found in ( 12T|) . The key property is that the denominator of the integrand consists of 
factors of (^X\Qi\X , where no two Qi are the same, along with some power of the factor 

{X\K\X . 

We find it helpful to rearrange the integrand in order to identify the cuts of 
master integrals as given in ( 123|) . ( 1271) . and ( 13T|) . This task is accomplished by partial 
fraction identities that split the denominator factors and reduce the power of (x\K\X 
if necessary. In effect, it is a reduction technique for the cut integrals. 

The splitting of factors with partial fractions proceeds as follows. First, split the 



factors 



among themselves, with the following identity: 



n.tlK- A] 


k 

= E- 

i=i 


1 




aj\Qi\X^ 


nil (Aig.iA] 


(A|g,|A 


1 im=l,my^i 


{X\QM~X) 



(41) 



Next, reduce the power of (fli^lf] in the remaining denominators, since the master 
cuts contain at most one: 





l[aj X] 


u;=l 


a,|g|A) 
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(^X\K\X 


"(A|g|A 


{X\KQ\Xr-' (^x\K\X 


(A|g|A 



(42) 



n-2 

E 

p=0 



n"=f ^ KI<5|A)) K-p-i 



m[n::Lp[at ~x] 



X\K\X 



p+2 



{X\KQ\Xy 



n—p—1 



Power-counting arguments ensure that enough appearances of A in the numerator to 
implement these identities as often as necessary. 
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It remains to implement the residue theorem, with the help of a generalized version 
of the differentiation identity flM|) and the careful treatment of higher-multiplicity poles 
arising in the factor {X\KQ\X). 

The procedure can be performed in generality. Formulas for the coefficients are 
given below. 

3. 3. Solutions for coefficients 

We now list formulas for the coefficients of master integrals [321 |33l [Ml [351 [36], through 
0{e). These references include D-dimensional versions of these formulas for going 
to higher orders in e, and with possible scalar masses. These generalizations will be 
described in sections to follow, but the more general formulas will not be reproduced 
here. 

Our starting point is a general form for a single term (i.e., the numerator is a 
monomial in the loop momentum) in the sum of the cuts of one-loop Feynman integrals, 

^(4vr)^ / (0^^''^(^') ^^^^((^ - (43) 

Here N is defined as the degree of T^^\i): it is the degree of the numerator minus 
the degree of the denominator, after setting £^ to zero everywhere it appears. In 
renormalizable gauge, N < 2. For convenience we also define some notation for the 
propagators, 

Q, ^ - K,. (44) 



\\Ki\X K^\\Qi\\ 
Di{e) = {i- Kf) = - K^A J = \ J . (45) 



\\K\\ {\\K\\ 

where the last equations are valid within the cut integral. Finally, the cut integral 
should be recast in terms of the spinors A, A with the replacement 

I = tx~x, 

and the value of t taken from (1^ . 



X\K\X 

3.3.1. Box coefficients The coefficient of the box identified by the two cut propagators 
along with Dr and Ds is given by 

C[Kr, K,, K] = \ {r^''\E)Drms{^))l_,p^^^;,^p^^^ + {Psr,l ^ Psr,2} (46) 

Here we use the following definitions. The vectors P^r.i and Psr,2 are the null linear 
combinations of Qr and Qg. 



Psr,l = Qs+{ ' ] Qr. (47) 
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Qs 



-Qs ■ Qr — V^s 
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{Qs ■ Qrf - QlQl- 
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(48) 
(49) 



3.3.2. Triangle coefficients If < — 1, the triangle coefficients are zero. If > — 1, 
the coefficient of the triangle identified by the two cut propagators along with Dg is 
given by 

1 



C[K,,K]^ 



2{N + l)\^r\Ps^Ps,2f^' 

/ / 1 N+l 

X ±^ ( r(^)(£)L>,(£)/A|X|A 



(50) 



\-^Qs\,\^Ps,1-tPs,2 



r~>0 



Here we use the following definitions. The vectors Ps,i and Ps,2 are null linear 
combinations of Qs and K. 



A, ^{Qs-Kf -QlK\ 



(51) 

(52) 
(53) 



The effect of the multiple derivative in the parameter r, evaluated at r = 0, is simply 
to pick out a term in the series expansion. 



3.3.3. Bubble coefficient There is just one bubble in the cut channel K. li N < 0, the 
coefficient is zero. If > 0, the coefficient is 



q=0 \ r=l a=q 



{r;a-q;2). 
N,N-a > 



,(54) 



s=0 



where 
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■N 
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Figure 3. A quadruple cut puts four propagators on shell. It is a trivial integral 
isolating a single box coefficient. 



-bW 



6!(m + l)v/A;""^ {Pr,i 

{Mv\Pr,2 



X 



dr^ 



I /\ n .A\\h 



_1 Af+1 
l,\\QrVi\\)' i\\K\\ 



\ 



v 



n+1 



T=0 



Here ?7, r/' are arbitrary spinors; they should be generic in the sense that they do not 
coincide with any spinors from massless external legs. 



4. Generalized unitarity: unphysical cut channels 

Although one-loop amplitudes are constructible by their branch cuts in physical 
momentum channels, and we have seen how this construction can be carried out, it can 
be useful to solve for the coefficients of master integrals by studying other singularities. 

Unitarity cuts can be "generalized" in the sense of putting a different number 
of propagators on shell. This operation selects different kinds of singularities of the 
amplitude; they are not physical momentum channels like ordinary cuts and do not have 
an interpretation relating to the unitarity of the S-matrix. Here, it becomes essential to 
work with complexified momenta. 



J^..^. Quadruple cuts 

The most direct application of generalized unitarity is to use a "quadruple cut" to find 
any box coefficient [37]. If we cut /owr propagators — equivalent to specifying a partition 
{Ki, K2, K^, K4) of the external momenta — then the four- dimensional integral becomes 
trivial. See Figure [3l 

^^^l-loop _ j ^4^ ^(^2) ^(^2) ^(^2) ^(^2) ^tree^tree^tree^tree ^55^ 

Applied to the master integrals, the quadruple cut picks up a contribution from 
exactly one box integral, namely the one with momenta {Ki, K2, K^, K4) at the corners. 
Therefore, the cut expansion collapses to a single term: 

^^^i-ioop ^ ^^^^^^ K,)AJ^{K,, K2, K,, K,). (56) 
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The quadruple cut of the scalar box integral is a Jacobian factor which is equal on both 
sides of the equation. The result for the coefficient is simply 

C4 = -J2 ''(£)Af<'(£), (57) 

where S is the solution set for the four delta functions of the cut propagators, 

5 = {£|£2 = o, i£-K,y = 0, {(i-K^-K^f = Q, (£ + 7^4)2 = 0}. (58) 

There are exactly two solutions, provided that momenta are allowed to take complex 
values. This is the origin of the factor of 2 in the denominator of (1571) . Thus it is easy 
to get all the box coefficients. 

In effect, the substitutions in the formula given previously in equation fH6l) 
implement the quadruple cut solutions directly. 



Triple cuts 

A triple cut, in which three selected propagators are put on shell, targets the coefficient 
of the unique scalar triangle including those three propagators. However, it also picks 
up contributions from all scalar boxes with those same three propagators. Moreover, 
there is a one- dimensional integral left over, so the calculation is less direct than the 
quadruple cut. 

Forde's parametrization of loop momentum allows the extraction of triangle 
coefficients from triple cuts [SS], for the case of massless propagators. Suppose the 
triple-cut conditions are given by 

£2 = 0, (£-Ki)2 = 0, {l + K^f = Q. (59) 

From Ki and i^s, construct null vectors as follows: 

where Si = Kf, = K^, and 

5'3(5'i-7) , 5'i(S'3-7) 
(y. = 

S1S3 - 72 ' S1S3 - 72 

Then the loop momentum can be expressed in terms of a single parameter t such that 
it satisfies the three cut conditions fl59l) explicitly. The parametrization is 

i = Kl + K',+ ^- {K'r\Y\Kt) + {KlnriKl'-) . (62) 

Plugging this expression into the integrand executes the triple cut, giving us the triangle 
of interest along with some boxes. Each box contribution has an additional propagator, 
with two poles in t. Therefore the box terms can be removed by expanding the expression 
around the point t = oo. The result for the triangle coefficient is given by 

- [lnUAiA2A,]{t)\^^,, (63) 



^ = Ki.K,±,/iKi.K,y-KfKl a = ^P^^, ^ = ^f^^. (61) 
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averaged over the two solutions for 7, where Ai,A2,A^ are the three tree-level 
contributions analogous to the amplitudes in f lTSj) . and the function Inf denotes the 
series expansion around t = 00. 

Another approach to triple cuts |l39j is to write the third delta function as the 
difference of two propagators with different ie prescriptions, 

2m5{p') ^ (64) 

and treat the problem as a double cut in the first two delta functions with modified 
propagators appearing in the tree-level amplitudes. Propagator modification is also 
the key to Feynman's Tree Theorem, relating loop amplitudes to sums of products of 
tree amplitudes produced by on-shell delta function insertions. The concept fits in 
the framework of generalized unitarity, so we will describe it briefly in the following 
subsection. 



4-3. Feynman's Tree Theorem 



Feynman observed [lUl HI] that any (Feynman) diagram with closed loops can be 
expressed in terms of tree diagrams, in the same spirit as Cutkosky but with more 
general cuts. It is not exactly the type of "generalized unitarity" described above for 
triple and quadruple cuts, because the propagators have been changed. We illustrate it 
in the simplest case, with scalar propagators at one loop. 

The key identity relates Feynman propagators Gp and advanced propagators Ga, 
which differ by a delta function. (Similarly, one could choose to use retarded instead of 
advanced propagators throughout.) Let k = {ko,k) be the 4-momentum of the scalar, 

m its mass, and u = \J {kl"^ + nn?- Then 

GF{k) 

GAik) 

Using the identity 
1 



'\k\ 
i 

i 

2w 



1 



kQ 



UJ + ie 
1 



k^ + u — ie 
1 



^2 



UJ 



ie kr. + UJ 



le 



uj^ + le 
i 



^2 



UJ^ 



ie sgn(A;o 



(65) 
(66) 

(67) 



where PV denotes the principal value prescription, we can see that 

GA{k) = GF{k)-'2i^6^^\k^ -m^). (68) 

Now consider the loop integral in which every propagator is replaced by an advanced 
propagator. All the poles are located above the real axis of w, so we can close the 
contour below, and the result is zero. But when we apply the identity fl68l) . we get 

0= /"Ar(fc) JJg5^(A;) = I N{k)^{GF{k)-2Txb'^^\k^ 

J k ■ J k 



m 



'■)) ■ (69) 



Here A^(A;) denotes the numerator. Expanding the product on the right-hand side, the 
first term is the original physical loop integral, and every other term has at least one 
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on-shell delta function. Feynman's Tree Theorem is this expression of the physical loop 
integral in terms of tree diagrams derived from its on-shell cuts. 

Recently, this theorem has been been used in the proof of covariance of the MHV 
prescription at one loop |12]. It is interesting to try to compute amplitudes directly from 
the tree-level ingredients with altered propagator prescriptions |l3l HU HS]. Compared 
to other "unitarity" methods, those based on Feynman's Tree Theorem work without 
reference to master integrals, as the entire integral is contained in the proper description 
of the propagators and their analytic continuations. 

4.4- Single cuts 

Quadruple cuts compute box coefficients. Triple cuts compute triangle coefficients, and 
box coefficients as well. Traditional double cuts compute bubble, triangle, and box 
coefficients, as we saw in the previous section. Single cuts should then give all the 
coefficients of master integrals. 

It is not easy to compute single cuts analytically. They are relatively unphysical 
and more strongly divergent than double (or triple or quadruple) cuts. When 
a single propagator is cut, the integrand is typically a tree- level amplitude with 
singular kinematics. Exceptionally, the kinematics are nonsingular for massless 
particles in supersymmetric gauge theories and some massive particles with additional 
supersymmetry [33], and they have been used to study higher- loop integrands in Af = 4 
SYM [16]. 

Nevertheless, some QCD amplitudes have been computed from single cuts alone 
[37] : single cut integrals have been analyzed and a prescription for their evaluation 
formulated in analogy with the residue theorem for double cuts [H]; and a single-cut 
method has been found using the principles of Feynman's Tree Theorem but grouping 
the terms together into a smaller set of single-cut amplitudes [44] . 

5. D-dimensional unitarity 

We have seen how to evaluate four-dimensional cuts to solve for the scalar box, 
triangle, and bubble coefficients in the expansion of a one-loop amplitude. In a four- 
dimensional expansion, there are additional cut-free "rational" terms. Traditional 
reduction techniques start with a loop integral and reduce it to a set of master integrals 
using linear relations generating the rational coefficients. None of the master integrals 
is purely rational. Where, then, do the purely rational terms come from? The answer 
is that they arise in the expansion in e. Higher-order corrections to the coefficients 
combine with the UV-divergent parts of the integrals to produce a purely rational term 
at C(eO). 

^i-ioop^^ E c„(K)L^o/.(K) + 7^ + 0(e) (70) 

n K={Ki,...,K„} 
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If we keep full dependence in e without the expansion, there are no terms apart from 
the master integrals in the original expansion 

^i-ioop^^ E c„(K)(e) X 4(K). (71) 

n K={Ki,...,K„} 

This property is the basis of D- dimensional cut-constructibility [HI [6l O [50]. We need 
to generalize our cut integral to D dimensions to pick up the higher-order corrections 
to the coefficients. Here we present the physical double-cut integral generalizing the 
treatment in Section |3] [511 152] . A similar analysis for generalized cuts is available as 
well [53] . 

All external momenta remain four-dimensional. We decompose the loop momentum 
i into its four-dimensional part i, and the remaining /j, which is orthogonal [54] 
ETJEO]. Then the integration measure factorizes, 



d^-''i = J d-''fi J dH = J du u-'-' J dH. (72) 

We have introduced the dimensionless variable m, defined by 

" = |r- (73) 

In the D-dimensional unitarity approach, we proceed by working on the familiar 
four-dimensional integral, which now depends on the new variable u. We keep the 
M-dependence throughout. Dropping the prefactor {At^Y /V{—e), the cut measure is then 

' du u~^~' [ d^l6^+\I^ - /i2)5(+)((£ - Kf - fi^) (74) 



To make use of the spinor technology that is so effective for four-dimensional cuts, we 
want to relate £ to a null 4- vector. This can be done with the help of the cut vector K. 
We let 

I=p + zK (75) 

where p should be null. This condition is implemented by a delta function on p'^ together 
with an integral over z. 



du u-'-' / dz / dH 6^+\f - fi^)6^+\{i - Kf - fi^)6{p^) (76) 



Now we have three delta functions. After some manipulation of their arguments, we are 
able to perform the 2;-integral trivially, finding that 



2 

The remaining cut measure is 



' - (77) 



^ du u~^~' j d^p S{p^) S{VT^ K^-2K-p). (78) 
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At this point we can replace the null momentum with spinor variables as before, with 
p = tXX. The cuts of the master integrals are modified as follows from the four- 
dimensional versions in Section |3l 

1 /•! , 

A/2 = / du u-^-'Vl-u (79) 

2Tn Jq 



A/3 = 

with 



2m 



du u 



-l-e 



:ln 



V^3 



u 



(80) 
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Kj ■ K 



B 



Aij, 




(83) 



Finally, we add the scalar pentagon to the basis of master integrals. Here we have three 
momenta /^j, Kj, K^ appearing in the three uncut propagators, and we define Qi, Qj, Qk 
by 



Qr 

The cut is given by 

A/5 = 



- {VT^)Kr + 



'l-Vl^){Kr-K) 



K^ 



K. 



du u 



{K^y 



i5) 



(86) 



S[Q3, Q2, Qi, K] ^^Q,.Q,- ^{Q, . Q,y - QIQ 



WiQ3 ■ - QlQl Q3-Q2 + ViQs ■ Q2] 



Q3Q2 



, SjQs, Qi, Q2, K] Qs-Qi- y^iQs ■ QiY - QlQl 



^^{Q3-QiY-QlQ\ • Qi + ^{Q, ■ QiY - QlQl 
S[Q2, Qi, Q3, K] Q2 ■ Qi - ViQ2-Qiy-QlQl 



WiQ2-Qiy-QlQl Q2-Qi + ViQ2 ■ Qi? - QlQl 
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where 



f K-Qk Qi-K Qj-k\ 



S[Qi, Qj, Qk, K] 



2det Q,-Qk Q- Qi-Qj 

\^ Qj ■ Qk Qi ■ Qj Q'j j 
^ Qk Qi' Qk Qj ■ Qk ^ 



(87) 



det Qi ■ Qk Qi Qi ■ Qj 
\^ Qj ■ Qk Qi ■ Qj Q"] I 



The full cut amplitude produces similar functional forms. As in the four- 
dimensional case, the pieces coming from different master integrals are identifiable 
by their unique logarithms. The pentagons are an exception: the arguments of their 
logarithms are the same as from the boxes. The cut pentagons are distinguished from 
cut boxes just by the rational prefactor in the functional form S'[Qi, Qj, Qk-, K\. In our 
approach, we use the logarithmic function to extract the box and pentagon coefficients 
simultaneously as total "box" coefficients, and later separate the pentagon parts from 
the box parts. 

The cut amplitude is naturally decomposed into linear combinations of the functions 
above, the cuts of master integrals, with coefficients rational in the momentum 
invariants — and now with additional dependence on the parameter u. This dependence 
turns out to be purely polynomial [M]. This behavior is expected on physical grounds 
but is not immediately apparent due to the irrational factor a/1 — u now appearing in 
the formulas. 

Because the w-dependence of the coefficients is polynomial, we do not need to 
finish the u integral explicitly. Instead, we use recursion/reduction formulas, derived 
from integration by parts, to deal with the higher powers of u. More specifically, terms 
in the cut amplitude will emerge naturally as combinations of the following functions. 








These functions are related to the cut master integrals by 



Bub^") = Fg^Bub^o) 

lliP\Z) = FS3(Z)Tri(°)(Z) + FS2(Z)Bub(°) 

Box(") = Fi^^BoxW + {Fi53(Zi)Tri(°)(Zi) + ^^(^^(^OBub^") + (Z, ^ Z^)} 
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{n-e).ZA 2^^[aJ {k-l/2-e). > 

Here {x)n = r(x + n)/r(x), and j = 2,3. In the equation for the box shift, Zi, Z2 
correspond to the two possible cut-triangles obtained by pinching the uncut propagators 
of the box. 

Alternatively, the result can simply be given in terms of dimensionally shifted 
master integrals. 

6. Treatment of massive particles 

In the preceding sections, we have explored unitarity cuts of amplitudes whose internal 
lines are all massless. When internal lines can come from massive particles, the formulas 
need to be modified. 

Scalar masses are fairly straightforward to deal with after having studied D- 
dimensional cuts. In dimensional regularization, the quantity /i^ behaves as a scalar 
mass for the four-dimensional parts of the loop momenta. It is not difficult to generalize 
the formulas for cuts and coefficients to account for different values of masses [32l [58| [35] . 
We can immediately calculate the coefficients of bubbles, triangles, boxes and pentagons. 

A bigger issue is that new types of master integrals arise: the scalar tadpole, and 
"massless bubbles," which have a single massless external leg on one side (K^ = 0) and 
a massive particle circulating in the loop. 

The scalar tadpole integral is 

/, ~ (92) 

e(e- 1) 

2 

777 

= _+rn\l-j- log(m2)) + 0(e), (93) 
and the massless bubble integral is 

^^^^ = °^~ Mf-M| e{l-e) ^^^^ 



; + (1 - 7) - MTZrjpiMf \og{Mt) - \og{Mi)) + 0{e) (95) 
1 ' ■ -r2^' 



or - - (7 + hgMf) + 0{e) if Mi = M2. (96) 

In massless theories, both of these integrals are set to zero in dimensional 
regularization. Neither of these integrals has a branch cut in any physical channel. 

Another problematic integral is a bubble with a single massive fermion on one side, 
which is also cut-free in four dimensions, proportional to (m^)"*^. However, the diagrams 
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contributing to the master integral expansion contain integrals that do have physical 
cuts in the channel, as well as contributions from self-energy insertions on an external 
line. One proposed resolution involves applying a unitary method in a specific gauge so 
that the wavefunction renormalization factor takes a fixed value [59]. 

The coefficients of the cut-free integrals can, in some cases, be completely fixed by 
the universal ultraviolet and infrared divergent behavior |60l EH [501 Ell ES] . 

Within the unitarity method, there is a proposal to find massless bubble coefficients 
from taking a suitable limit and the tadpole coefficients from an artificial double cut, 
constructed by introducing an auxiliary, unphysical propagator into the formulas [M] . 
based on the integrand classification of [65]. The artificial double cut can be evaluated 
by standard unitarity methods. 

Finally, within generalized unitarity, techniques for computing single cuts 
analytically can be used to derive some QCD amplitudes [13 HE]. Single cuts exhibit 
some pathological behavior in general and need to be treated carefully. The subject 
deserves further study. 

7. MHV diagrams 

A discussion of new techniques for loop amplitudes would not be complete without 
mentioning the applications of MHV diagrams by Cachazo-Svrcek-Witten (CSW) rules 
as a tool in amplitude construction. 

MHV diagrams and CSW rules [30] are an alternative to and on-shell analog of 
Feynman diagrams and rules. The external legs are in fixed helicity states. The vertices 
are MHV amplitudes, and there is a prescribed on-shell continuation for propagators 
to reconcile momentum conservation and intermediate states of definite helicity. At 
tree level, it is understood that they are valid because they reproduce all the correct 
singularities. 

The MHV diagram technique can be applied in several ways at loop level. Most 
immediately, it can be used to construct tree-level amplitudes for use in a (possibly 
generalized) unitarity cut integral, and the four-dimensional cuts evaluated as above. 

Similarly, cut propagators that have been analytically continued to lie on shell can 
be continued further to D = A — 2e dimensions; one can then continue with explicit 
evaluation of a D-dimensional cut integral, which was used in [661 E3 EHl E9] to recover 
one-loop MHV amplitudes in Yang-Mills theory with various amounts of supersymmetry. 
Extensions of CSW rules to include massive scalars [70] lead to MHV-diagrammatic 
constructions of rational terms in one-loop gluon amplitudes [71]. 

The CSW rules were conjectured based on the localization of Yang-Mills amplitudes 
in twistor space. Along these lines, twistor space actions were written that reproduce 
the MHV diagram construction [721 [73]. The twistor actions are compatible with, but 

lack, the h amplitude at tree level and the all-plus amplitude at one-loop. These 

amplitudes are not constructible from MHV diagrams but are also well localized in 
twistor space [ZH [30] . 
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Transformations have been found converting the standard Yang-Mills or QCD 
Lagrangians into Lagrangians whose tree-level expansion reproduces the CSW rules 
[75| 176] . relying on an expansion in lightcone gauge. One would expect that the 
transformed Lagrangian then includes the exceptional interaction vertices (— — + 
amplitude at tree level and the all-plus amplitude at one-loop) in some form. Indeed, 
they show up as counterterms after suitable regulators are introduced [771 [ZH] • However, 
the regulators either break Lorentz invariance or violate the equivalence theorem and 
exclude supersymmetry. 

8. Recursive constructions of coefficients or rational terms 

BCFW-type on-shell recursion relations can be constructed for parts of loop amplitudes 
lacking branch cuts. Most obviously, the rational terms of one-loop amplitudes seem 
suitable for such treatment. It is also interesting to consider recursion relations 
constructed for the coefficients of the master integrals, as these are also purely rational 
functions. 

Indeed, if we calculate the coefficients from cut integrals, then they are naturally 
related to the tree amplitudes on either side of the cut. Recursion relations for these tree 
amplitudes lead to recursion relations for the coefficients. The difference is of course that 
we now consider shifts involving the loop momentum, and therefore the pole structure 
is somewhat more subtle. 

Recursion relations are derived from factorization properties. Study of coUinear and 
multiparticle factorization limits shows that it is safe to apply a BCFW shift on two legs 
attached to the same vertex of the master integral, as long as this shift is valid within 
that particular tree amplitude [H]. With such a shift, kinematic poles involving the 
loop momentum remain unaffected. Recursion relations of this type have been written 
explicitly for gluon amplitudes in a split helicity configuration, with either an A/" = 1 
chiral multiplet or a complex scalar circulating in the loop. 

On-shell recursion relations can be used to compute the rational parts of massless 
one-loop amplitudes [791 [80l [HI [82] (see [83l[8l] for reviews). We can get the coefficients 
of master integrals purely from four- dimensional cuts, and find the remaining parts with 
the residue theorem directly in four- dimensional kinematics as well. 'Rational term" 
could be an ambiguous name, since the master integrals contain some rational terms 
along with their logarithms; here we define it as the remaining four- dimensional part of 
the amplitude when all scalar integrals are set to zero, as in (iTOl) : 

^^^i-ioop_^ ^ c„(K)L^o^.(K) (97) 

n K={Ki,...,K„} 

We assume that the coefficients of the integrals have been computed first. 

Now, examine the poles of this rational part. They include not only the expected 
physical poles, but also spurious poles that will be cancelled in the sum with the cut part. 
An example of this kind of unphysical singularity can be seen in the denominator of the 
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two-mass triangle; one can verify that the singular behavior is rational. These spurious 
poles can be simple poles or double poles. Moreover, it is sometimes difficult to find a 
momentum shift such that the amplitude vanishes in the limit that the shift parameter 
z goes to infinity. Recursion relations for rational parts of one-loop amplitudes have to 
deal with both of these problems. 

The spurious poles can be canceled explicitly by redefining the rational term. In 
equation ( 1971) above, we have, schematically for an m-point amplitude, 

A^"'°°P = C„ + 7^^ (98) 

where Cm is the cut-containing part with master integrals and TZm is the rational 
term. There is a systematic (but not unique) way to "complete" the cut part by 
adding some rational terms so that the spurious poles are cancelled, starting with two- 
mass triangle integrals and including higher- dimensional box and triangle integrals as 
necessary [851186]. If these extra terms are represented by CTZm, then the cut completion 
is 

Cm Cm CT^m (99) 

and we recast the rational part as 

nm = nm- CTZm. (lOO) 

Now a recursion relation can be constructed for TZm from knowledge of its physical poles. 
In a given factorization limit, there are three types of contributions, corresponding to 
the location of the original loop integral. The recursion relation takes the form 

Yl ('^LiZa)^AR{z^) + ^^.(z^) + Al{zS-^Ar{zM ■ (101) 

poles a ^ " " a / 

Here, Al and Ar are lower-point tree amplitudes; TZl and TZr are rational parts of 
lower-point one-loop amplitudes; and the "factorization function" J-" is derived from 
diagrams where the loop is a correction to the propagator that is going on shell [85] . 

This three-term recursion relation is still not completely correct, because there 
are generally also boundary terms coming from the non-vanishing of the amplitude as 
z — 7- oo. And when it is possible to select a shift such that the amplitude does vanish at 
infinity, then there is a similar issue in the recursion relation itself: along with the three 
recursive terms in equation f llOip . there are extra terms in the residues themselves that 
do not fit the recursive pattern. 

The resolution is to combine information from two different shifts. First, use a 
primary shift, under which the residues are purely recursive as in equation (IIOII) . Then 
look at the boundary term, described by the large-2; behavior of the amplitude and 
denoted by Inf ^4;'^^°°^. This boundary term can be constructed from a recursion relation 
using an auxiliary shift. The auxiliary shift may have non- recursive terms in its residues, 
but they can be dealt with as long as they vanish in the large z limit. We also need 
to compute Inf Cm, the large z behavior of the cut-completed terms, under both shifts 
separately, and take care to avoid double- counting the residues at physical poles in the 
cut-completed and rational parts. 
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9. Survey of recently published one-loop amplitudes 

Since the introduction of the unitarity method in 1994, many new analytic, explicit 
expressions for one-loop amplitudes have become available. Here we list those that are 
the most complete (all helicity configurations and rational terms). 

Gluon amplitudes in = 4 SYM are known explicitly for MHV [IT] and NMHV 
[87] configurations, which is enough to give all results through seven external legs [88] . 

In QCD, partial results for five-parton amplitudes [HHl |90] were completed with the 
unitarity method and summarized in [91] 

The amplitudes for e~^e~ — )■ four partons (also necessary for W/Z/Drell-Yan 
production in association with two jets at hadron colhders, or three-jet production 
in deep inelastic scattering) have been given in [92], which also featured the used of a 
triple cut to find a three- mass triangle coefficient. 

The components of the six-gluon amplitude [56l|9a[I7l[Tl|92IMll28l|80l[T8l[29l 
|85llHIll95] have been summarized in [961 197] with an accompanying Mathematica code. 
The results have been cross-checked with a semi-numerical derivation [98] . 

The six-photon amplitude in QED has been found from both multi-cut and form 
factor decomposition methods and cross-checked in [99] . The result was given for scalar 
QED and Af = 1 QED by cut methods in [lOOj . Both analytic forms confirmed the 
numerical results of [101] . 

The components of the Higgs + four parton amplitudes have been found in the 
limit of large top quark mass [TDa [TDal HDl HDSl ^OQi [Ml UDE]. The analytic results 
give a fast (10 ms) evaluation of the full cross section, and the semi-numerical result 
[1091 1110] has been confirmed. 

Amplitudes featuring massive fermions are a bit more rare. Four-gluon scattering 
via a massive quark loop was fully calculated using cut methods in [50] . 

Formulas for the ttgg amplitude have been found from cut methods and described 
in [63 11111] and have been checked against the earlier (less compact) results of [11211113] . 

Lately, the calculation of the py-boson mediated process — )• duQQii [114] has been 
presented, an analytic version of [115] and including the spinor correlations present in 
W-boson decay. 

10. Beyond one loop 

Unitarity and recursive methods have been applied to compute multiloop amplitudes, 
but almost only in theories of maximal supersymmetry, which merit a separate 
discussion. 

Maximal supersymmetry implies many internal cancellations and hence a relatively 
simple integral structure. Moreover, Af = 4 SYM and JV = 8 supergravity possess 
symmetries even beyond those directly implied by gauge symmetry and supersymmetry, 
simplifying the results even further. 

It is worth noting that even in maximally supersymmetric theories, multiloop 
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calculations have been completed only for MHV configurations. No one should be 
astonished, then, that progress with reduced supersymmetry is taking more time. 

With on-shell methods, there are several obstacles to address to proceed from one 
loops to two or more. First, some convenient simplifications within QCD do not carry 
over. At one-loop, we could apply a color decomposition and calculate color-stripped 
"partial amplitudes." The subleading-color contributions were related by permutation 
identities to the leading-color part |116[ [T7] . At two loops, we can start with a color 
decomposition, but the subleading contributions need to be computed separately. 

In QCD, one-loop gluon amplitude calculations have routinely proceeded according 
to a supersymmetry decomposition: first, complete the internal gluon with a full A/" = 4 
supermultiplet, and compute the amplitude in the A/" = 4 theory; then, remove the 
new fermion contributions again by subtracting four times the A^ = 1 amplitude 
with a chiral multiplet in the loop; finally, cancel all the scalar by explicitly adding 
a nonsupersymmetric amplitude with a complex scalar in the loop. The decomposition 
is written schematically as 

AAf=0 gluon _ AjKf=4: A AAf=l chiral , ^scalar (TCtO^ 

^1-loop ~ ^1-loop ~ ^^1-loop "T ^1-loop- l-LU^J 

The underlying principle is that supersymmetric amplitudes are simpler due to internal 
cancellations (for example, diagrams loop amplitudes in TV = 4 SYM have no internal 
triangle subdiagrams), and that a nonsupersymmetric amplitude with a complex scalar 
loop is simpler than with a gluon loop, because there are no spin degrees of freedom. At 
two loops, there is no relation of this kind. One can hope to get some insight by starting 
with a calculation with internal scalars, which can help in formulating an ansatz for the 
full solution [117] . but it is not a systematic approach. 

Even more seriously, in general theories, the master integrals are more numerous 
and less explicit; and a unitarity cut of an amplitude is unavoidably D-dimensional, 
typically with internal spins. The four-dimensional spinor-helicity formalism has been 
tremendously useful in writing amplitudes analytically. Higher-dimensional spinor 
formalisms have been explored |118l 11191 1120[ I12H I122j but still need further attention. 
Continuations to larger integer values of D are enough to fix all rational terms of one- 
loop amplitudes |123] and certainly help in defining the tree-level input for cut-based 
methods. 

Two-loop master integrals can be generated for any specific process by the Laporta 
algorithm |124[ 1125] , which is based on integration- by-parts identities relating various 
integrals in systems of (large numbers of) linear equations. Programs implementing 
the Laporta algorithm efficiently are available |126l 11271 1128] . Recently a framework 
has been outhne for reducing two-loop integrals to a finite basis avoiding doubled 
propagators, whose singularities are more difficult to control [129] . It still remains 
to list them explicitly, check independence, and evaluate them. 

The operations that have been proposed to construct one-loop amplitudes from 
their single cuts using different propagator prescriptions |l3l HI] can be extended to 
higher-loop integrals |15], since there is no need for a basis of master integrals, and 



Loop amplitudes in gauge theories: modern analytic approaches 



29 



the information regarding integration can be said to be carried in the propagator 
prescriptions. 

In QCD, unitarity methods have been apphed to compute the 4-gluon amphtude 
to two- loop order |117] , as well as the two- loop splitting amplitudes (universal functions 
governing collinear behavior) for g ^ gg [130] and q ^ qg and g ^ qq |131] . Seven 
master integrals were needed, of which four vanished through 0{e) and two others 
were products of one-loop integrals. In principle there could have been two-loop planar 
integrals that had no cuts in two-particle channels; fortunately, these did did not appear 
in this particular calculation, and so it sufficed to look at cuts in the physical channels. 
The corresponding amplitude with internal scalars were found first as an exercise. 
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